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ABSTRACT 

Planetary rovers provide a considerable challenge for robotic systems in that they must operate tor 
long periods autonomously, or with relatively little intervention. To achieve this, they need to have 
on-board fault detection and diagnosis capabilities tn order to determine the actual state ot the vehicle, 
and decide what actions are safe to perform. Traditional model-based diagnosis techniques are not 
suitable for rovers due to the tight coupling between the vehicle's performance and us environment. 
Hvbrid diagnosis using particle filters is presented as an alternative, and its strengths and weaknesses 
are examined. We also present some extensions to particle filters that are designed to make them 
more suitable for use in diagnosis problems. 


1 INTRODUCTION 

Planetary rovers provide a considerable challenge for robotic systems in that they must operate tor long periods 
autonomously, or with relatively little intervention. To achieve this, they need to have on-board fault detection and 
diagnosis capabilities in order to determine the actual state of the vehicle, and decide what actions are safe to pertorm. 

The diagnosis problem is to determine the current state of a system given a stream of observations of that system. In 
traditional model-based diagnosis systems such as the Livingston system [l] that was a part o t e emo ^ . = 
experiment [2,3] diagnosis is performed by maintaining a set of candidate hypotheses (in the Remote Agen exp . 

SETS* hictois was Lpfl abou, fe cufTcni ««e of .ho systtfn. mi using the mod.Uo pr.d.a lh. «pb»d 
future state of the system given each candidate. The predicted states are then compared with the observations of what 
actually occurred. If the observations are consistent with a particular state that is predicted, that state is kept as a 
candidate hypothesis. If they are inconsistent, the candidate is discarded. Traditional diagnosis systems typica y use ; a 
logic-based representation, and use monitors to translate continuous-valued sensor readings into tscr ' 
variables. The system can then reason about the discrete variables, and compare them with the predictions of the model 
using constraint propagation techniques. 

Unlike spacecraft, rover performance depends significantly on environmental interactions. The on-board sensors 
provide streams of continuous valued data that varies due to noise, but also due to the interaction between the rover and 
its environment. For example, a rover may have a sensor that reports the current drawn by a wheel. In normal opermmn. 
this quantity may vary considerably, increasing when the vehicle is climbing a hill, and decreasing on ownw s pc . 
t£ diugnoL systeni feetis to be able to tiistioguish t, change in the onren, Jruwu due to the termtn being traversed 
from a change doe to a fault in the wheel. A second issue for rovers is that their weight and power is very tightly 
constrained. For this reason, any on-board diagnosis system must be computationally efficient, and should^ able to 
adapt to variations in processor availability. Ideally, we would also like it to adapt based on its own performance, 
spending more time on diagnosis when a fault is likely to have occurred, and less time when the system appears to be 

operating normally. 

A rover's dose coupling with its environment poses a considerable problem for traditional model-based diagnosis 
systems. A particular sensor reading may be normal under certain environmental conditions, but • ndicat. v e o f a fault 

others, so anv monitor that translates the sensor reading into a discrete value such as nominal , or o -no * 

must be sophisticated enough to take all the environmental conditions into account. This can mean that the ^agnems 
problem is effectively passed off to the monitors— the diagnosis system is very simple, but re . 

values from extremely complex monitors that 'diagnose" the interaction between the system and us nv.ronment as p 
of translating continuous sensor values into discrete variables. To overcome this problem, we need to reason directly 
with the continuous values we receive from sensors. That is. our model needs to be a hybnd svs.e.n. consis ing of a s« 
of discrete modes that the svstem can be in. along with a set of continuous state variables. The dvnam.es ot system 


is do'wnbed m lerrns >1 a set * »l ei|tialii>MN governing E he evolution ot the »tutc varumes. nut these eqnainnis will fie 
■ liit'ereni m different modes. In addition. a transition function describes -how the o, tern moves Irom one mode to 
. u M 1 1 her. and an observation function defines the likelihood ot an observation given 'tie mode and die values ol the 
ovsiem variables. 

Plus Inhrid model is verv similar to a partittllv t>hsrr. able \Uirkt>v dec: sum proevw (POMDPh POM DPs are Ireijuently 
used as a representation lor decision-theoretic planning problems, where the task is ’<> determine the best action to 
perform given the current estimate of the actual state of the system. This estimate, referred to as the relief Mute, is 
exactly what we would like to determine in the diagnosis problem, and the problem of keeping the beliet ^tate updated 
is well understood in the decision theory literature. The belief state is a probability distribution over the system states, 
and is updated by “pushing it through" the system model to produce a new probability distribution over the possible 
future stale of the system, and then conditioning that distribution on the actual sensor values that were observed to 
remove anv states inconsistent with the observations, reduce the probability of any state from which the observations 
are unlikely, and increase the probability of states that predict the observations well. We will take essentially this 
approach to the diagnosis problem in this paper. 

Unfortunately, the algorithm we described above for maintaining an accurate beliet state is computationally intractable 
for the types of problem we are interested in. Since our model contains both discrete and continuous variables, the beliet 
state is a set of multidimensional probability distributions over the continuous state variables, with one such distribution 
for each mode of the system. This distribution may not even be unimodai, so just representing it is a complex problem, 
but applying the model to it to predict future state is infeasible due to the severe computational restrictions on-board a 
rover. Therefore, an approximation needs to be made. The approach we will take is to use particle filtering [4,5]. 

A particle filter represents a probability distribution using a set of discrete samples, referred to as particles , each of 
which has an associated weight. The set of weighted particles constitutes an approximation to the belief state, and has 
the advantage over many approximation methods that it can represent arbitrary distributions. To perform diagnosis, we 
apply the predictive model to each particle individually, and then condition on the observations by multiplying the 
particle's weight by the probability of observing the sensor readings if the particle represented the true state of the 
system. To prevent a small number of panicles from dominating the probability distribution, the particles are then 
resampled, with a new set of panicles, each of weight one, being constructed by selecting samples randomly based on 
their weight from the old set. We will discuss panicle filters in more detail in Section 3 below. 

Particle filters have already proven very successful for a number of tasks, including robot navigation [6], Unfortunately, 
they are less well suited to diagnosis tasks. This is because the mode transitions that we are most interested in 
detecting — namely transitions to fault states — are those with the lowest probability ol actually occurring. Thus, there is 
a risk that there will be no particle in a fault state when a fault occurs, and so the system will be unable to diagnose the 
fault as no particle will predict the observations well. To apply particle filters successfully to rover diagnosis, we will 
have to overcome this problem, and we discuss two approaches here. Tn the first, we increase the number ot particles 
when the current particles seem to be doing a poor job of predicting the observations. The idea is that when the total 
weight of the particles falls, it indicates that no particle is predicting well. In response to this, the diagnosis system can 
suggest that a fault may have occurred, and use additional resources (additional particles) to try to identity the tault. 

A more promising approach is to think of the particle filter as a convenient way to divide the available computation 
time between individual states. In this model, the system uses as many particles as computational limitations will allow, 
and divides them up between the system modes according to how well each mode predicts the sensor readings as 
before, but also according to how important it is to diagnose each mode. So tor example, it a particular mode m predicts 
the observations well at present, we would ensure that there are always some particles in a tault state m that is 
reachable from m with non-zero probability. This means that il the fault did occur, there are guaranteed to be some 
particles in tn' that would predict the observed behaviour after the fault. 

In the next section we discuss the rover model in detail. In Section 3 we describe particle tiltenng and demonstrate its 
weaknesses when applied to diagnosis problems. This paper describes a work in progress, so in Section 4 we will 
describe our proposed modifications to the standard particle filter in detail, and in Section 5, present some preliminary 
results on real rover data, using a simple version of our proposed approach. The final section looks at the relationship 
between this work and some previous approaches to this problem, and discusses some future directions tor this work. 


modklmm; a planktary roykr 


\ ; , wo .aid above. wo mode! a rover as a hybrid .vstem, Nio discrete oompononl ol (ho rovers >tato ropiosonls the 
various operational and fault modes of the rover, while l he continuous >late describes the speed ot the wheels, the 
current heme drawn hv various subsystems, and mj on. Following | |. we model a iover as a tuple <\l, V. T,E.O> where 
the elements of the tuple are as follows: 

* XI is the set of discrete modes the system can he in. We assume that V/ is finite, and write m tor an individual 


system mode. 

• V' is the set of variables describing the continuous state ot the system. 

• r is a transition function that defines how the system moves from one mode to another over time. We write 
Pr iim.m') for the probability that the system moves from mode m to mode m . We may also include a second 
transition function Pr, im.ajn’) which is used when an action ,i occurs. This gives the probability of moving 
from m to tn when act ton a is executed. 

• E is a set of equations that describe the evolution of the continuous variables over time. The equations that 
apply at a given time potentially depend on the system mode, so we write E m tor the equations that apply in 
mode m. These equations will in general include a noise term to account tor random variations in the *tate 
variables. Here we will assume Gaussian noise, with the parameters ot the Gaussian determined individually 
for each equation. 

• 0 is a function mapping the system state into observations. We will assume that the observable system 
characteristics are some subset of the system variables V, with their values corrupted by Gaussian noise (again 
with parameters that may be a function of the variable, and the system mode), so we write 0(v,m) for the 
observed value of some variable v in mode m. 

[n addition to these, we will write Pr(s'\s) for the probability distribution over future states s’ given some state s\ where 
y and s are hybrid states, so Prfs'js; includes both the distribution over the future mode given by Pr?/™ \m) % and the 
distributions over the continuous variables given by E m . 


The diagnosis problem now becomes the task of determining the current mode m that the system is in, and the values of 
all the state variables in V. We will assume for the rest of this paper that alt the state variables are observable, that the 
problem of tracking their values through the noise is relatively straightforward (if the noise in the observations 
overwhelms the true values of the variables, we have little chance to determine the system mode), and assume that the 
diagnosis problem is to determine the mode only. Note that the particle filter approach we describe does determine the 
most likely value for the system variables as well as estimating the mode. 

For a system as complex as a rover, we would in general take a component-based approach to modelling, where each 
component (for example, a wheel) is modelled separately and the complete model is the cross-product of its 
components. Since each component has much simpler dynamics, and typically a smaller set ot discrete modes ot 
operation, this makes the modelling task much easier, but it also simplifies diagnosis by reducing the number of faults 
that need to be considered to only those in components that aren’t behaving nominally. Of course, we can expect a 
certain amount of dependency in a component-based model (for example, we would expect each wheel to be moving at 
approximately the same speed in most circumstances). Taking advantage ot these dependencies allows us to simplify 
the diagnosis problem in some cases, for example by identifying a wheel that is faulty because its behaviour doesn t 
match that of the other wheels. However, we will not discuss these issues in detail here. 


The experiments we will present in Sections 3 and 5 use actual telemetry data from NASA Ames’ Vlarsokhod rover 
The Vlarsokhod is a planetary rover built on a Russian chassis that has been used in field tests from 1993-99 in Russia, 
Hawaii, and the deserts of Arizona and California. The rover has six independently driven wheels, and for the 
experiments we present here, the right rear wheel had a broken gear, and so rolls passively. The Vlarsokhod has a 
number of sensors, including ones that measure body and arm geometry, and battery current. We will restrict our 
attention to diagnosing the state of the broken wheel, and will therefore use only data from the wheel current and wheel 
odometry sensors. We will treat each wheel independently in the diagnosis. For each wheel, we have a model, taken 
from [7], with the following characteristics (in fact, the model of each wheel is identical): ^ 

• ’ Xt consists of 23 system modes, including for example an idle state, the wheel current rising after a command 

is given, uniform speed driving, and a total of 14 different fault states. 

• V consists of variables for the wheel current and wheel speed, and the derivatives ot current and speed. 

• T is a fairly sparse matrix, with at most six successors for any given mode. The probability of a transition to a 

fault state is 0.01 or less. We include a separate transition function for the start and end of a command. 



( Y(-:itc .i ,et >1 n particles vliere each particle p has i .tale n and a w<:, -hl A. ts 
from the initial state di arihutnui. ami /\ =/. 


dimpled randomly 


I -or :ach lime step. do: 

a. For each particle p. do: 

i. Select a future state p./ for p hv sampling from Pr ( ^ l p ), the disirihution over possible 
future states given by the model. 

ii. Re-wemhl p by main plying Us weight by the probabi iiiv i )t the observations <) eiven p, as 
follows: 


b. 


p\ - Pr ( o | p*)p* 

Resample ;i new particles p, by copying the current particles where each particle p is added to the 
new samples with the following probability: 

Pn p. = p) = ~- r 

yp- 


Figure (: The particle tillering algorithm. 


• The state equations in E consist of the previous value plus a constant term and noise. The noise is Gaussian 
with standard deviation in the range 0.001 to 1.0. and the equations are independent for each state variable. 

• The equations in O are independent for each variable (but vary depending on the mode), and consist only of 
the variable’s value plus Gaussian noise. Again, the standard deviation varies trom 0.001 to 1.0. 

3 PARTICLE FILTERS 

A particle filter approximates an unknown probability distribution using a weighted set of samples. Each sample or 
particle consists of a value for every state variable, so it describes one possible complete state the system might be in. 
As observations are made, the transition function is applied to each panicle individually, moving it stochastically to a 
new state based on its current dynamics, and then the observations are used to re-weight each panicle to reflect the 
likelihood of the observation given the panicle's new state. In this way, panicles that predict the observed system 
performance are highly weighted, indicating that they are in likely states of the system [6]. The major advantage of 
particle filters is that their computational requirements depend only on the number of panicles, not on the size of the 
system. This is of huge importance to us as it allows us to do diagnosis in an anytime fashion; increasing the number of 
panicles when there is computation time available, or when we suspect a fault has occurred, and decreasing the number 
when other operations require the processor. To implement a panicle filter, we require three things. 

• A probability distribution over the initial state of the system. 

• A model of the system that can be used to predict, given the current state according to an individual panicle, a 
possible future state of that panicle. Since T is stochastic, and E includes noise terms, the predictive model 
selects a new state for the panicle in a Monte Carlo [8] fashion, choosing by sampling from the probability 
distribution over possible future states. 

• A way to compute the likelihood of observ ing particular sensor values given a state. In our case, this is given 
by the observation function O. 

The panicle filtering algorithm is given in Figure l. Step (i) is the predictive step, where a new state is calculated in a 
Monte Carlo way for each panicle, and this new state is then conditioned on the observations in step (ii). In the case of 
our rover model, the initial distribution is uniformly in the Idle mode, in which there is no power to any of the wheels. 
The predictive step is performed by applying T to each panicle, and then applying the appropriate equations from £ to 
the state variables, sampling values from the Gaussian error terms. Once the panicles have been re-weighted, we can 
then calculate the probability of each mode simply by summing the weights ot the panicles in the mode. 

Problems with particle filters for diagnosis 

Unfortunately there are a number of difficulties in applying particle filters to diagnosis problems. In particular, the 
filter must have a particle in a particular state before the probab.lity of that state can be evaluated. It a state has no 
particles in it. the assumption is that its probability of being the true state of the system is zero. This is a problem in 
diagnosis problems because the transition probabilities to fault states are typically very low. so panicles are unlikely to 
end up in fault stales during the Monte Carlo predictive step. Without a particle in a fault state, that state will continue 



to he zr-cn /.cru probability ot hjtnm ■K.airreil, 
known is .1 itnplc trnptnert \hrrn‘nt. 


il no other aate is dome a iiood job ot predicting 


!'hts filiation is 


Figure 2 illustrates 'his problem. Each graph allows the most likely modes that the *heel is in. shown over part ol one 
,,f the trials m which the .v heel is initially idle, and then at .tep 12 is commanded to drive forward at a Used speed. The 
araphs on the left .how the performance of Wheel I. which is operating nominally. Hie graphs on the right show the 
performance of Wheel 6. which is faulty. In the top line, the probability of the fault occurring is 0. 1 rather than its true 
value oft). t)|. Here the fault is quickly detected til Wheel 6. The middle line ot graphs diows the same situation with the 
fault probability set to 0.05. Here there is a significant delay between the wheel being commanded, and the detection ot 
the fault. In the bottom line of graphs, the fault probability is set to its true value, and in this case the tault is not 
successfully detected because insufficient panicles enter the fault state. One might expect that once a particle enters the 
fault state its weight would be high since it would predict well, and at the re-sampling step it should lead to several new 
particles being created. Unfortunately, this did not occur in this situation because although some particles did enter J-he 
fault state, their continuous parameter values did not agree with the observations well, so they still had low weights. The 
continuous parameters did not match because each of the particles that entered the fault state came from the 
CmmandedRunnin g state, in which the current and wheel speed are expected to be much higher than the observed 
values. In the next section we examine two techniques for overcoming this problem and ensuring that faults are detected 
even when [heir probability ot occurring is low. 


4 USING DOMAIN KNOWLEDGE TO DIRECT THE PARTICLE FILTER 

The simplest solution to the sample impoverishment problem is to increase the number ot particles being used. Given 
the constraints imposed on on-board systems, this approach is probably unrealistic. The data presented above used 
10.000 particles per wheel, and runs in Java in approximately 0.5s per update on a 750MHz Pentium 3. This is probably 
at the upper limit of the number of panicles we could expect to use on-board a rover — the time available for diagnosis is 
longer, but the computation will be much slower. Thus running with ten times as many particles (which is rough y 
equivalent to multiplying the fault probability by ten) is probably impractical on the rover, and even 10.000 particles 
may be unrealistic as the model gets more complex. This could be somewhat overcome by only increasing the number 
of particles when there is some evidence that the system is predicting poorly. In order to achieve this, we need some 
measure of when this occurs. The obvious measure is to look at the total weight of the particles after conditioning on the 
observations. If no particles are predicting the observations well the total weight should drop. Unfortunately, in practice 
this is rarely useful because there are a number of other possible causes for this behaviour. For example, 
moving from a state in which there is high confidence in the sensor readings to a state with more sensor noise will tend 
to drop in weight even if they are still predicting the observations well. We see this in the Marsokhod model because the 
idle mode has' relatively large variance for the observation noise, whereas the CommandedRunnmg mode has smal er 
variance, so the total particle weight increases when the system moves from the Idle to the CommandedRurming mode, 
even for wheel 6 where CommandedRurming predicts the observations poorly. 

Another way to reduce the likelihood of sample impoverishment is to bias the distribution ot particles in favour of states 
which correspond to important faults. In this approach, the particle filter no longer approximates the true be let 
distribution, but instead uses some oracle to identify •'important” future states, and allocates some portion of t e 
particles to those states. The particle filter can then be thought of as a way of dynamically allocating our available 
computation. Each system mode is assigned a fraction of the available computation time, in the form ot a number ot 
particles. As in the standard particle filter, particles move from state to state via Tand E, and via the resampling step, 
but the oracle ensures that every state it considers important is assigned at least some small traction ot the available 
computation. This ensures that if a transition to one of the important states occurs, there will be particles in the state so 
its probability will be non-zero. If the particles in one of these states predict the observations well (which they should it 
it is the true state), their weight will then increase, at which point the normal operation of the particle tiller increases the 
number of particles in the state. 

The question that remains is how to implement the oracle. For a complex system such as a planetary rover, with many 
components each with its own set of possible failure modes, there are exponentially many possible failure modes, so 
this is a non-trivial problem. However, one approach that seems promising is to use a traditional model-based diagnosis 
system such as Livingston [ I ]. Livingston operates much more quickly than hybrid approaches because it does not 
consider the continuous dynamics of the system. Instead, it uses a discrete abstraction ot the problem, and tries to 
identify likely states that are consistent with the observations. We pointed out in the introduction that this approach is 
not in general suitable for diagnosing rovers, but it is suitable for the purpose of identifying likely system modes to 
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Figure 2: Panicle filter results for wheels I (nominal performance) and 6 (faulty, with a broken gear). In the top row 
the probability of the fault is ten times its true value, in the middle row it is 5 times, and the bottom row has the true 
fault probability. The fault is easily detected in the top row, detected alter a delay in the middle row, and not detected 

in the bottom. 


apply some computation in. By asking Livingston to provide a number ot candidate hypotheses that explain the 
observations, we can with high probability be confident that the true mode is in one ot the hypotheses. We then assign 
some computation to all of the modes consistent with the hypotheses. The speed ot the qualitative system should mean 
that using the oracle would not significantly impact computation, which still should depend only on the number ot 
particles. The integration of Livingston with the particle filter approach is currently work in progress, as it adds a 
number of additional complications including building an additional system model, and ensuring that the discrete and 
hybrid models agree with one another and can easily be translated back and torth. 
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Figure 3: Results for particle filter with bias. All states with > 25% probability were used as starting points tor the 
forward search, and 0.5% of the particles were assigned to each of the found states. On the left are results tor wheel 

l, and on the right for wheel 6. 

For simpler systems such as the Marsokhod wheel diagnosis we have used in this paper, the Livingston- based approach 
is unnecessary. Instead, we can use an oracle based on forward search from the current high-probability states. Since 
each system mode in this model has at most six possible successors, and there are typically only two to three high 
probability modes at any time (often one of these succeeds another), we find in practice that in most cases a simple one- 
step look-ahead search adds fewer than five modes to those that already contain particles. 

5 RESULTS 

The results we present here are based on the Marsokhod model we described in Section 2. Dr. Rich Washington 
supplied the model and the data, which came from his work on using Kalman filters tor rover diagnosis [7]. The only 
changes made to the model were to make it suitable for use with a particle filter; no changes were made to model 
parameters or transition probabilities. To demonstrate our approach we use a small piece ot one of the telemetry data 
files (the same piece used in Section 3) in which the rover is initially idle, and then a drive command is issued, resulting 
in an increase in current to each wheel, followed by a corresponding increase in speed, and then a constant speed. As 
before, wheel 6 is faulty, with a broken gear (this corresponds to the Gear&EncoderFaultRunning state in the model). 


Figure 3 shows the results for the biased panicle filter. For these results we used single step forward search from all 
states with probability > 0.25 to select the set of bias states. Each of these states was then guaranteed to receive at least 
0.5% of the total number of particles at each re-sampling step. The left hand graph is the probable states tor Wheel 1, as 
before. Like the sraph in the bottom row of Figure 2, the Premature Action state was given high probability before step 
13. This state appears in situations where the effects of an action are seen before the signal to perform the action is seen, 
due to problems with the rover telemetry. In this case it is a spurious result due to the model of the Idle state not 
allowing sufficient noise in the observations. A small adjustment to the model would remove this problem, which is 
only present in the data for two of the wheels. The right hand graph shows the same data tor Wheel 6. In this case, the 
fault state is found at step 20. seven steps after the command to drive the wheel was observed. This compares 
favourably with the results in Figure 2, where the fault was detected three steps after the command in the top row (ten 
times probability), and 10 steps after it in the middle row (five times probability). 

6 DISCUSSION AND RELATION TO OTHER WORK 

One closely related piece of work is Verma et al’s decision-theoretic particle tiller (9|. The problem they are attempting 
to solve is essentially the same one we are interested in. but their approach is to assign a utility to every state and use the 
utility to change the distribution they draw the particles from. Effectively, (his is equivalent to altering the transition 
function st) (hat the probability of a transition from ->tate s with utility m.vi to state v with utility u{ * ) becomes 


I’rr. . ittr.L This has a amilar effect m mcreasum ‘he probabilities as we did in Retire 2. l-nra-laliveU Minplo 

ill amu >si . i asks uich as the one we have presented here, the approaches seem /erv similar. However, lor more complex 
i asks, we believe that tisme a discrete diagnosis tool as an oracle in direct the computation ol the particle biter will 
allow us to make more effective use ol the available computation than the utility based method, which will increase the 
probability ol transitions to a potentially very high number of fault stales, especially as any reasonable utility function 
would give all fault slates a high utility. 


Another related effort is the work of Washington [7] that applies Kalman Filters to this problem. In this work, the 
continuous dynamics in each mode is tracked by a set of Kalman filters. The mam problem with this approach is that 
the number of (liters tends to increase over lime because each time a transition is made to a state the initial conditions 
for the filter are different, and filters with different initial conditions cannot be combined. This is not a problem for 


particle filter-based approaches because the particle (liters can represent arbitrary distributions over the parameter 
values. >o particles entering a state with two different sets of initial conditions will form a bi-modal distribution. As we 
said above, we used the model and data from this paper in our own work. We see fewer errors in the mode identification 
with our approach than in Washington’s paper, although we are sometimes slower to identity the fault, and our 
computational requirements are somewhat higher. 


Lemer et al [10] use linear-Gaussian Bayesian networks to represent a hybrid system, and perform diagnosis by 
computing the true belief state at every step. This avoids the problems we have described, which are directly caused by 
the approximation approach. However, the computational requirements are far too high to be used on-board a rover. 


As we said in the introduction, this is a work in progress. There is still much work to do on the problem of how to 
integrate a model from Livingston with this system to act as an oracle. We have demonstrated that a simple look-ahead 
search performs quite well, but this is clearly inadequate for large diagnosis problems. We are also examining a number 
of other approaches to improving diagnosis with particle filters, such as backtracking when prediction is poor, and re- 
sampling past states based on observations that occurred more recently. Finally, we are investigating how a diagnosis 
system of this type would fit with the CLARAty rover architecture [111 being used for future NASA missions to Mars. 
We are actively involved in developing these approaches to be integrated with the Mars’07 mission. 


We would like to thank Dr. Richard Washington for providing the model of the Marsokhod rover we used tor this 
paper, along with the telemetry files from his field tests with the rover. This work has also been influenced by 
discussions with Dr Washington, Vandi Verma, and Dr. Robin Morris. 
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